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I"""- We consider the experimental implications of the down type isosinglet quark, D, predicted by the £5, group 

to Higgs boson searches at the LHC. The pair production of D quarks at the LHC and their subsequent decays 
D — > hd and D — » Wu has been analyzed. For a light Higgs boson of mass 0(120 GeV), an analysis based on 
fast simulation of the ATLAS detector response shows that the bb channel becomes as efficient as the yy channel 
^ , for discovering the Higgs particle if mo < 630 GeV. 
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I. INTRODUCTION 



^vq ■ Although the standard model (SM) explains the results of the particle physics experiments performed up to now 

J> I with an accuracy high enough to withstand even the challenges set by the precision measurements, it still leaves 

some of the basic questions open. The number of elementary particles, the reason for their mass hierarchy, and 
the unification of gravity with the other known forces are some examples of these open issues. The grand unified 
theories (GUTs) aim to answer at least some of these questions by imposing a fundamental symmetry between all 
known fermions of the same family. This symmetry, manifesting itself at high energies, is expected to reduce the 
, number of free parameters in the SM. The experimental implication of extending the existing SU"c(3) x SUw(2) X 

■ Uy(l) group structure of the SM to a single gauge group with a large fundamental representation is the prediction 

of new particles. The exceptional Lie group £g has been long considered as one of the favorite candidates for 
such a GUT gauge symmetry group 01 0]. The new colored particles predicted by £g are isosinglet quarks, 
leptoquarks and diquarks depending on model variations. Some of these particles may soon be accessible by the 
LHC experiments |3j,|4|]. The current limit on the mass of the down type isosing let quark is m D > 199GeV (.5]. The 
upgraded Tevatron could reach mo ~ 300 GeV Jfl, whereas the LHC will cover the region up to mo ~ 1 TeV for the 
pair production channel [7] independent of the mixing between d and D quarks. For the single production channel, 
if the sine of the mixing angle exceeds 0.025, the discovery reach would be mo ~ 1.5 TeV 18Q. In this note, we 
consider the possible impact of the mixing between the isosinglet quark and the SM down type quarks on the 
light Higgs boson h (e.g. to/, < 135 GeV) searches with the ATLAS detector yfl. According to the ATLAS technical 
design report |3], for such a light Higgs boson, the h — ► yy channel is the only one that would give a possible 
discovery with more than 5 (J significance with 100 fb _1 of integrated luminosity. However, if the model is true 
and the mass of the lightest additional quark is suitable, its decays involving the Higgs boson enhance the Higgs 
discovery potential. In particular the channel studied in this note, pp — ► DD —> hjW j —> bbjlvj, becomes as 
efficient as the yy channel for discovering the Higgs particle if mo < 630 GeV. 

This paper is organized as follows: Section [TTJ contains the description of the model and the new interactions 
involving the Higgs particle. Section|III]briefly discusses the outcomes of the new interaction: A) the possibility of 
differentiation between the quarks of this model and any other model with an additional down type quark and B) 
an increase in the Higgs production rate at the LHC. Section QV] elaborates on the second outcome by describing 
the performed Monte Carlo (MC) studies based on a fast simulation of the ATLAS detector response. Section IVl 
summarizes the statistical analysis that estimates the experimental reach to discover a light Higgs for the early days 
of the LHC. 



II. THE HIGGS INTERACTION 



If the SM SV~c(3) x St%(2) x E/y(l) group structure originates from the breaking of the £g GUT scale down to 
the electroweak scale, then the extended quark sector should be written as: 
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As shown, each SM family is extended by the addition of an isosinglet quark, respectively denoted by the letters D, 
S, and B. The mixing between the new and the SM down type quarks is responsible for the decays of the former. 
In this study, the intra-family mixing of the new quarks is assumed to be dominant compared to their inter-family 
mixing. Considering a similar mass hierarchy between the SM and the new quarks, the D (B) quark is assumed to 
be the lightest (heaviest) one. The Lagrangian relevant for the D quark weak interactions is : 



L D = 



2\/2sin 9w 
4sin0w 



u 6 Ta (1 - y 5 ) d cos (p + u e Y a (1 - y 5 )Dsm<j) 



W a 



(2) 



sin cos (j) 
cos 9w 



dy a {\-Y 5 )D 



\ / 47ta e , 



12 cos 9w sin 6\y 



[Dy a (4 sin 2 9 W - 3 sin 2 <j> (1 - y 5 )) D + dy a (4sin 2 0^-3 cos 2 0(1 - 75)) d] Z a + h.c. 



where <p is the mixing angle between the d and D quarks, and 9 represents the usual CKM quark mixings, taken 
to be in the up sector for simplicity of calculation. Assuming that the squared sum of the row elements of the new 
3x4 CKM matrix gives unity, the measured values of V ud y us , and V u t, JH constrain sin0 < 0.07 0]. The total 
decay width and the contribution from neutral and charged currents were already estimated in As reported in 
that work, the D quark decays through a W boson with a branching ratio of 67% and through a Z boson with a 
branching ratio of 33%. 

The origin of the masses of SM particles is explained by using the Higgs mechanism. The Higgs mechanism 
can also be preserved in group structure as an effective theory, although other alternatives such as dynamical 
symmetry breaking are also proposed ifioll . On the other hand, the origin of the mass of the new quarks (D, S, B) 
should be due to another mechanism since these are isosinglets. Therefore, the mass terms in the Lagrangian, 
written in the mass basis, are simply: 
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where rru and mo are the masses of the d and D quarks respectively. Based on this initial consideration, the 
studies in ||7l|9|l did not consider the channels involving the Higgs boson h for the D quark decays. However, the 
mixing between d and D quarks will lead to decays of the latter involving h after spontaneous symmetry breaking 
(SSB). To find these decay channels, the interaction between the Higgs field and both down type quarks of the first 
family should be considered before SSB. In this study, the Lagrangian in SM basis, containing also the cross terms 
between d and D with mass-like coefficients (a f /o, m D( i), is written as: 
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The four interaction coefficients (a dd , a d jj, mjj d , muo) in equation © can be obtained in terms of the two quark 
masses and the two mixing angles between the left and right components of the two down type quarks. In order to 
use the new variable set, the first step is to write the mixing between the SM and the mass states as follows: 
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where <j) L / R is the d-D quark mixing angle with the subscript L (R) standing for the left (right) components of the 
quark fields and the superscript (M) standing for the SM (mass) basis. After the SSB, the Higgs field is expanded 
around the new minimum, v, as H = v + h to separate the terms that contribute to interaction and mass coefficients: 
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Figure 1: The decay width of the D quark for sin0 = 0.045: The lower curve corresponds to the case without Higgs interaction, 
whereas the upper curve has been calculated by taking into account the Higgs channels for mu = 120 GeV . 
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where v = tj / \/2 and 77 = 247 GeV is the vacuum expectation value of the Higgs field. After simple calculations, 
the Lagrangian for the interaction between d, D quarks, and the Higgs boson becomes : 



L m = '^W(p L D M D M h 
sin0 L cos0 L ^ M 



(7) 



(l-7 5 )m D + (l + 7 5 );7v 
(l + 7 5 )m D + (l-7 5 )m f/ 



d M h 



D M h 



2v 

sinfccosfe - M 
2v 

+ ^± ^ <p L d M d M h 

V 

Note that the right mixing angle has completely disappeared from the final formula 0. The total width of the 
D quark as a function of its mass is shown in Fig. [T]for the illustrative value of sintp = 0.045. It is seen that the 
D quark has a narrow width and becomes even narrower with decreasing values of since it scales through a 
sin 2 <p dependence. The relative branching ratios for the decay of the D quark depend on both the D quark and 
the Higgs mass values. For example, at the values of D quark mass around 200 GeV and the Higgs mass around 
120 GeV: Br(D -» Wm)~60%, Br(D -> W)~12%, Br(D -> ZJ)~28%, whereas as the D quark mass increases the 
same ratios asymptotically reach 50%, 25% and 25% respectively. As the Higgs mass increases from 120 GeV, 
these limit values are reached at higher D quark masses. The branching ratios as a function of the D quark mass 
are given in Fig. |2]for two values of the Higgs mass: 120 GeV (solid lines), which is just above the limit imposed 
by LEP II results, and 135 GeV (dashed lines), where the bb decay mode is taken over by the W + W~mode. 
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Figure 2: The D quark branching ratios as a function of the D quark mass for two Higgs mass values: 120 GeV & 135 GeV. 

Table I: For pair production of D quarks, the decay channels involving at least one Higgs boson. The branching ratios and the 
number of expected Higgs particles are calculated assuming m//=120 GeV and m£>=250 (500) GeV. 
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III. THE EXPERIMENTAL IMPLICATIONS 



Studies about the D quark pair production without the Higgs particle have been previously reported elsewhere 
0]. This section concentrates on the possible impact of the interactions involving the Higgs particle. For example, 
equation ^} shows that, the Higgs boson could be used in the s channel for the creation of D — d quark pairs: /; — ► 
Dd or dD. However, as the Higgs couplings depend on the quark masses, the right handed D quarks (Dr) would 
be produced {'^) 2 times more often than the left handed ones (Di). Due to the left handedness of the charged 
weak interactions, the produced D quarks would subsequently decay mostly via the neutral current channel. The 
precise measurement of such a process could help in differentiating the £g model from other models that also 
involve additional down type quarks. 

Depending on the masses of the D quark and the Higgs boson itself, the model could also boost the overall 
Higgs production at the LHC. This boost is particularly interesting for the Higgs hunt, one of the main goals of 
the ATLAS experiment. For example, if the D quark mass is as low as 250 GeV, its pair production cross section 
at the LHC becomes as high as 10 5 fb as can be seen in Fig. [3] Convolution of the cross section with the D quark 
decay branching ratios from previous section (17% for D — > hj at itid=250 GeV) yields the number of expected 
Higgs events. In the low mass range considered in this note (from 1 15 up to 135 GeV), the branching ratio h — * bb 
is about 70% JH]. The signatures of the expected final states involving at least one Higgs boson are listed in table 
|I]together with the yearly expected number of events for two example D quark masses: 250 and 500 GeV. For the 
remaining of this note, although the case involving the Z is more suitable from the event reconstruction point of 
view, we will concentrate on the last row, which has the highest number of expected Higgs events per year. 
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Figure 3: Tree level pair production of D quarks at LHC as a function of the new quark's mass. 



IV. MONTE CARLO STUDY 



The Lagrangian presented in section |TT] has been implemented in a tree level event generator, Comphep 
v4.4.3.|l 1;1, to investigate the possibility of detecting the Higgs particle and reconstructing it from b— jets us- 
ing the ATLAS detector. Assuming a light Higgs boson of mass 120 GeV, four mass values for the D quark have 
been taken as examples: 250 GeV, 500 GeV, 750 GeV, and 1000 GeV. 10000 signal events were produced for 
each mass value under study with the W h j j (h — > bb &W — > £v) final states using the CTEQ6L1 PDF set lfl2ll . 
Other parameterizations, such as CTEQ5M and CTEQ5L1 were also tested to investigate the systematic effect of 
the PDF selection on the DD pair production total cross section. It was found that the difference was less than 
10% over the range of 200 to 1600 GeV, CTEQ6L1 giving the lowest results up to a D quark mass of 1 TeV The 
generator level cuts on the partons, guided by the performance of the ATLAS detector, are listed as: 

K\ < 3.2 , 
p Tp > 15GeV , 
R p > 0.4 

where rj p is the pseudo-rapidity for the partons giving rise to jets; pjp is the transverse momentum of the par- 
tons; and R p is the angular separation between the partons. The imposed maximum value of r\ p requires the jets 
seeded by the partons to be in the extended barrel region of the calorimeter where the jet energy resolution is 
optimal. The imposed lower value of pj p ensures that no jets that would eventually go undetected along the beam 
pipe are generated at all. The imposed lower value of R p provides good separation between the two jets in the 
final state. Using the interface provided by CPYTH v2.3 IU3I1 . the generated particles are processed within the 
ATLAS software framework version 11.0.41. The decay and hadronization processes were handled in PYTHIA 
ifTill and detector response simulation was performed with fast Monte Carlo program ATLFAST [15 ] which uses 
a general parameterization of the detector, ignoring the individual event topology. The default properties of the 
parameterized fast simulation are discussed elsewhere Q[l5l]. One should also note that in this work, the recon- 
structed b— jet energy and momenta were re-calibrated like in ^ to have a good match between the mean value of 
the reconstructed Higgs mass and its parton level value. 

As for the background study, all the SM interactions giving the W jb jt j j final state have been computed in 
another tree level generator, MadGraph v2. 1 . |16il, using the same parton level cuts and PDF. The SM background 
cross section is calculated to be 520 ±1 lpb. The reasons for using two separate event generators, their compatibil- 
ity, and their relative merits have been discussed elsewhere Qj, The generated 40000 background events were also 
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Figure 4: The distributions of kinematic parameters for SM background and signal for an example D quark mass of value of 
500 GeV. The solid histograms are for signal and the dashed ones for the background events. The leftmost plot contains the 
Pt distribution for the light jets which show clear distinction between signal and background cases. It also contains the pj 
distribution for for the i-jets (dotted histogram) which has the same shape for both signal and background. However, the middle 
plot shows that the fo-jets originating from the signal events can be discriminated from the background ones by considering the 
angle between the two jets. Finally the rightmost plot contains the scalar sum of all transverse momenta providing another 
criterion for separating the signal from the background. 



processed in the same way using the ATLAS software framework for hadronization and calculation of detector 
effects. 



A. Analysis Details for mo = 500 GeV 

The W boson will be reconstructed using an electron or a muon and the missing energy associated with the 
undetected neutrino, whereas the b— jets will be combined to reconstitute the h boson. These two bosons can be 
merged together with the remaining two jets to obtain the invariant masses of the two D quarks. Therefore, the 
final state particles of interest are j/, jb jj&Er^iss f° r which the distributions of the kinematic values are shown in 
Fig- El On the leftmost plot, one can observe that the pj distribution of the jets coming from D quark decays peaks 
around 200 GeV (solid line), whereas the same distribution for the background jets is much less energetic: the pj 
distribution of the b— jets peaks around 40 GeV (dotted line) and for other jets the peak is around 120 GeV (dashed 
line). This difference can be used to reject the background events most of which (about 80%) are originating from 
the tt production. In these events, one W decaying leptonically and the other hadronically give rise to the same 
final state particles as the signal. Another property that can be used to discriminate the background events is the 
fact that two b— jets in the signal originate from the same particle, whereas for most of the background cases, they 
originate from two different particles. The cosine of the angle between these jets, shown in Fig. |4]middle plot, is 
therefore peaked at one for the signal events, as opposed to the background case, where the distribution is more 
uniform. The scalar sum of all transverse momenta in an event (Hj), is given in the same figure, rightmost plot, 
with again a solid line for the signal events and with a dashed line for the background events. This variable will 
be of great benefit for distinguishing the signal events from the background, especially for the high D quark mass 
values. 

The reconstruction starts with the requirement of 1 electron or 1 muon and at least 4 jets, two of them identified 
as b— jets. For the leptons, the transverse momentum minimum value (pj.e > 15 GeV) is selected. The transverse 
momentum requirement on the two most energetic jets and the cosine of the angle between the two b—jet are 
optimized to get the best signal to background ratio for the final D quark invariant mass distribution. The h 
boson invariant mass is formed using the two b tagged jets remaining after the angular cut. To further reduce the 
background, the invariant mass of the two non b tagged jets is also calculated, and the events yielding a value 
smaller than Mz are rejected. As there are no Z bosons in this search channel, this cut ensures the purity of the 
selected events by removing the background events originating from the hadronic decays of the SM Z boson. 
Another rejection criteria applied at this stage is related to Hj, where only events with Hj > 800 GeV are selected. 
The association between the two most energetic jets and the reconstructed bosons is not unique. Both possibilities 
are calculated and the case with the smallest difference between the two reconstructed D quark invariant masses is 
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Table II: Optimized event selection cuts and their efficiencies for iri£> = 500 GeV. 
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Figure 5: Reconstructed invariant masses of the D quark (left) and of the Higgs boson (right) together with the SM background 
(dotted lines) after 30 fb _I integrated luminosity. The mass of the D quark is set to 500 GeV and Higgs boson to 120 GeV. 



selected. As the last consistency check, the difference between the two reconstructed D quark masses is required 
to be smaller than 100 GeV, a value compatible with the resolution of the reconstructed mass. The event selection 
cuts imposed on both signal and background are given in table [TT]in the order they were applied in the analysis. 
The same table also contains the individual efficiency values of each cut for both signal and background cases with 
an error of a few percent. As can be calculated from the table, the final efficiency for the signal is about 1 1 % and 
for the background is about 2%. 

The invariant mass distributions after all cuts for the D quark and the Higgs boson are presented in Fig. for 
30 fb _1 integrated luminosity. After defining the signal region for D as Mp ± 50 GeV and for h as M/, ± 30 
GeV, the number of events for the signal (S) and the background (B) can be integrated for both to obtain the 
statistical significance a = S/y/S+B. For this set of parameters, it is found that the D quark can be observed with 
a significance of 13.2(7 and at the same time the Higgs boson with a significance of about 9.5(7. One should note 
that, in the SM Higgs searches, such a high statistical significance can only be reached with more than 3 times 
more data: with about 100 fb _1 integrated luminosity. 



B. Extension to other mass values 



An analysis similar to the one presented in the previous section was performed for the other three D quark 
masses: 250, 750 and 1000 GeV. For each mass, the cut values were re-optimized to get the best statistical signif- 
icance in the Higgs boson search. The method is to start from a low value of the Hj cut and scan the available 
sample by increasing the cut value, like a sliding threshold to look for a peak in the reconstructed mass histogram. 
The invariant mass distribution of the reconstructed Higgs boson for different D quark mass values is given in Fig. [6] 
where the dashed line shows the SM background, the red data points are for the signal events only and the black 
data points are the sum of the signal and the background events. For each increasing D quark mass value under 
consideration, the total integrated luminosity was also increased to be able to observe the Higgs boson signal. 

Table HIT1 gives the number of events for signal and background processes after 3 years of low luminosity data 
taking time (30 fb integrated luminosity) for all considered D quark mass values. With such integrated luminos- 
ity, the 3(7 signal observation limit for the D quark is at 750 GeV. If mo > 1000 GeV, the required luminosities 
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Figure 6: The reconstructed masses for the Higgs particles for different D quark mass values. The dashed histogram represents 
the SM background, the red data points are for the signal case and the black data points are the sum of signal and the background. 
From left to right 250, 750 & 1000 GeV D quarks were considered using 10, 100 and 3000 fb -1 integrated luminosities 
respectively. 



Table III: The expected number of signal (S) and background (B) events for D quark and Higgs boson searches after 30 fb 
integrated luminosity. The counted events are in a range of ±50 GeV (±20 GeV) from the generator level value of the D quark 
(Higgs boson). 
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45.1 


27.6 
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9.5 


3.5 


2.1 


0.84 


0.42 



can be attained either by running the LHC for many years at the design parameters or by considering a possible 
upgrade, named as Super-LHC, which would yield 1000 fb _1 per year. On the other hand, the significance val- 
ues presented in Table iHTl could be reduced by systematic errors which have not been evaluated here, and by the 
dilution of the statistical significance which occurs in searches for an excess in any bin of distribution. 

V. CONCLUSIONS 

This work gives the Lagrangian for the D quark including its interaction with the known SM particles and Higgs 
boson. One of the two implications of the Higgs interaction is discussed in detail: the boost to the Higgs discovery 
potential at the LHC. It is shown that the ATLAS experiment can use the D quark decay channels to enhance the 
Higgs discovery potential independent of the mixing angle between d and D quarks. The prospects for the Higgs 
discovery are discussed using Monte Carlo and fast simulation techniques for a light Higgs particle (120 GeV) and 
some exemplary D quark mass values ranging from 250 to 1000 GeV. Fig. |7]contains the 3(7 (dashed lines) and 
the 5 <7 (solid lines) reaches of Higgs boson (triangles) and D quark (circles) searches for the above mentioned 
values. Therefore, a light Higgs boson could be discovered with a 5(7 statistical significance using the DD — > hW jj 
channel within the first year of low luminosity data taking (integrated luminosity of 10 fb ) if mo < 500 GeV. 
Under the same conditions but with one year of design luminosity (integrated luminosity of 100 fb -1 ), the 5 (7 
Higgs discovery can be reached if mo < 700 GeV. This is to be compared with the studies from the ATLAS TDR, 
where the most efficient channel to discover such a light Higgs is the h — > yy decay. This search yields about 
8(7 signal significance with 100 fb -1 integrated luminosity. The presently discussed model could give the same 
significance (or more) with the same integrated luminosity if mo < 630 GeV. Therefore, if the isosinglet quarks 
exist and their masses are suitable, they will provide a considerable improvement for the Higgs discovery potential. 
The same channel also provides the possibility to search for the lightest of the isosinglet quarks, providing a 5 (7 
discovery signal if mo < 800 GeV, within 100 fb _1 of integrated luminosity. However one should note that the 
version of the fast simulation software used in this note was not fully validated by the ATLAS collaboration. 
The study of other channels involving leptonic decays of the Z boson with thoroughly validated versions of the 
simulation and reconstruction software is in progress. 
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Figure 7: The 3 and 5 (7 reach of ATLAS in the Higgs and D quark search for increasing D quark mass values 
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